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We study theoretically a strongly type-II s-wave superconducting state of two-dimensional Dirac 
fermions in proximity to a ferromagnet having in-plane magnetization. It is shown that a magnetic 
domain wall can host a chain of equally spaced vortices in the superconducting order parameter, 
each of which binds a Majorana-fermion state. The overlap integral of neighboring Majorana states 
is sensitive to the position of the chemical potential of the Dirac fermions. Thermal transport and 
scanning tunneling microscopy experiments to probe the Majorana fermions are discussed. 



Majorana fermions are hypothetical particles that are 
their own antiparticles, as originally proposed to describe 
neutrinos [1 j. Because of the non-Abelian statistics as 
Ising anyons in two spatial dimensions [2] and their po- 
tential application to a fault-tolerant quantum computa- 
tion [2-4 , the Majorana fermions have attracted revived 
interest in condensed-matter physics [5]. Theoretically, 
they appear as zero-energy states bound to vortices in 
topological superconductors [6], including v = 5/2 frac- 
tional quantum-Hall states [2], s-wave superconductors 
of two-dimensional (2D) Dirac fermions [7], and chiral 
p-wave superconductors [6j [8|. However, it remains a 
challenge to realize and manipulate Majorana fermions 
in a practically detectable and controllable manner. 

Recently, there have been several proposals for re- 
alizing Majorana fermions as the edge modes at het- 
erostructure interfaces. It has been argued that the su- 
perconductivity induced by the proximity effect on the 
surface of a three-dimensional (3D) topological insula- 
tor (TI) [9HTT] can host Majorana modes propagating 
along the edge of the surface. When the superconduct- 
ing region is surrounded by two magnets having oppo- 
site out-of-plane magnetizations, two charge-neutral Ma- 
jorana edge modes are combined into a charged Dirac 
fermion, which can be probed with charge transport. A 
Rashba-semiconductor thin film sandwiched by an insu- 
lating magnet having an out-of-plane magnetization and 
an s-wave superconductor has also been proposed to pro- 
duce a chiral Majorana edge mode [T2HT5] . 

In this Letter we show theoretically that a magnetic 
domain wall in an insulating ferromagnet attached to a 
2D strongly type-II s-wave Dirac-fermion superconductor 
can produce a chain of vortices [Fig. [TJa)], each accom- 
modating a Majorana zero mode [3], when the magneti- 
zation points to the out-of-plane direction at the head- 
to-head magnetic domain wall. A hybridization of Ma- 
jorana fermions at neighboring vortices can be tuned by 
the chemical potential of the Dirac-fermion superconduc- 
tor, which can be probed with thermal transport and 
scanning tunneling microscopy (STM) experiments. We 
also propose two heterostructures that can realize such 



electronic states: (A) an interface of an insulating fer- 
romagnetic thin film and a bulk s-wave superconductor 
whose normal state is a 3D TI [Fig. [TJb)] and (B) a device 
structure with stacked s-wave superconductor, Rashba 
semiconductor, and insulating magnetic film [Fig. [IJc)]. 

We consider a 2D continuum BCS-Dirac Hamiltonian 
with a single Dirac cone in the single-particle dispersion, 

H = J d 2 xifj^ (x){[ivD x e 3 - B(x)] ■ a - /i<j }V ; (^) 

+ Jd 2 x{A(x)^(x)ia 2 ty\x)] T + K.c.}, (1) 

under the Zeeman field B(x), i.e., an exchange field 
coming from a ferromagnet. Here, v is the Dirac-fermion 
velocity, \i the chemical potential, and the covariant 
derivative D = V — ieA with A being the vector po- 
tential. The spinors ^{x) = (cL.,cL) and i/j(x) = 
( C £ct' c «|) crea ^ e an d annihilate fermions at position x 
with charge e, respectively. We use a = (01, 02, 03) for 
the Pauli matrices as well as 00 for the 2x2 unit ma- 
trix. The unit vector es is normal to the 2D plane, while 
x = xiei + #2^2 lies in the plane. In the absence of 
the magnetic field B(x) = 0, the superconducting order 
parameter A (as) is assumed to take a uniform constant 




FIG. 1. (Color online) (a) Schematic picture of the modulus 
of the superconducting order parameter (thin surface) in ei- 
ther of the two proposed heterostructures shown in (b) and 
(c). Along the boundary between two domains of a ferro- 
magnetic insulator (M) with opposite in-plane magnetization 
(bold arrows), a chain of vortices emerges. A single Majorana 
fermion is bound to each vortex (dots). 
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A , which we choose real and positive. Henceforth, we 
assume strong type-II superconductivity so that the or- 
bital coupling to the vector potential A can be ignored 
in the analysis of A(x). 

Here, we introduce a head-to- head domain wall that 
extends along the e2-direction [Fig. (2^a)]. Away from 
the domain wall, the Zeeman field is given by B — >• ±Bei 
for x\ — > ±00, respectively. Therefore, it can be approxi- 
mated with a uniform in-plane Zeeman field ±Be\. This 
shifts the Dirac cone and therefore causes a finite center- 
of-mass momentum q = ±ge2 of the spin-singlet Cooper 
pairs with q — 2B/v [16]. As a consequence, the su- 
perconducting order parameter A will acquire a spatial 
modulation A(sc) = A e iq ' x for x\ —> ±00. 

Alternatively, the spatial modulation of A(x) by 
the in-plane Zeeman field can be described using the 
Ginzburg-Landau (GL) free energy, 



F=Jd 2 x [a|A(^)| 2 + 7 |D A A(^)| 2 +^|A(x)| 4 ] 



(2) 



where V/\ = V — 2i(eA + v~ 1 A) and we introduced the 
real parameters a, 7 > 0, and u > 0. The linear coupling 
to the Zeeman field, written in terms of the fictitious 
vector potential A := es x B, is allowed due to the broken 
mirror symmetry about the 2D plane. Away from the 




FIG. 2. (Color online) (a) Magnetic moment of the ferro- 
magnet around the domain walls of types (i) and (ii). This 
structure continues in the other spatial direction, (b) and (c) 
Color maps of the phase of the superconducting order param- 
eter A(sc) in our numerical solutions of the GL model Eq. ^ 
on a grid of 50x50 sites for the types (i) and (ii) of head-to- 
head domain walls, respectively. They have been obtained at 
the criticality by ignoring the quartic term. In the e<i direc- 
tion, open boundary conditions and periodic boundary con- 
ditions were chosen for cases (i) and (ii), respectively. Fixed 
boundary conditions resembling the solutions deep inside each 
domain were imposed in the e\ direction. The domain- wall 
parameters were set to / = 6 lattice sites and s = 1/4 while 
B/v was chosen such that two periods of the phase modula- 
tion fit into the grid. A variation in the model parameters I, 
s and \B/v\ does not change the result qualitatively. 



domain wall, the mean-field solution in the strongly type- 
II limit can indeed be obtained as A(x) = A e zq ' x with 
A = ^/\a\/2u for a < 0. We employ this GL theory 
to determine a spatial profile of A (as) in the domain- wall 
region where the Zeeman field B varies with x. 

In fact, the profile of A (a?) strongly depends on 
whether the magnetic moment of the attached ferromag- 
netic film changes its direction along the domain wall (i) 
within the plane as in stripes of magnetic thin films [17] 
or (ii) towards the out-of-plane direction as in magnetic 
thin films with large magnetic anisotropy [18] . The mag- 
netic moment near these domain walls can be modeled 
as a divergence-free function, 



M(x,x 3 ) B 



\=±s 



e,\ — Ae 7 



tanh 



\x\ 



XI 



(3) 



with j = 2 for the case (i) and j = 3 for the case (ii), 
where I and s control the characteristic width of the do- 
main wall. We assume that the Zeeman field B acting on 
the Dirac-fermion superconductor at the interface x% = 
is given by B(x) = M(x,0). A key observation here is 
that the out-of-plane component Bs(x) vanishes. Then, 
this in-plane Zeeman field induces a fictitious magnetic 
field defined by B := V x A, which becomes B = and 
B = es(B/l) sech 2 (xi/l) for cases (i) and (ii), respec- 
tively. In analogy to the vortex-lattice phase of a type-II 
superconductor, we would thus expect to obtain vortices 
at the domain wall only for case (ii). This is confirmed 
by the numerical solutions of the phase of the supercon- 
ducting order parameter A (a?) for the cases (i) and (ii), 
presented in Figs. |2jb) and (c), respectively. For case (i), 
the contour of the phase of A (as) follows the magnetic 
flux B and hence no vortices emerge. For case (ii), a 
chain of vortices in A(x) appears along the domain wall. 
All the vortices have the same winding number +1 or — 1, 
where the sign is determined by the sign of B/v. 

The emergence of the vortices in case (ii) can also be 
understood from the integration of the phase of A(x) 
deep inside the two domains, i.e., along the contour 
drawn in Fig.[2|c). Provided A(x) does not vanish along 
the contour, the flux enclosed by this contour takes an 
integer multiple of the superconducting flux quantum $0, 
that monotonically changes by 1 each time the contour 
length along the domain wall is increased by ir/q. It 
indicates that the intervortex distance is given by ir/q. 
In Fig. [2|c) we observe that A (a?) changes sign when 
translating it by the separation of neighboring vortices 



x Q 



(0, ir/q), i.e., A{x) = -A(x + x q ), Va3. 



Having determined the spatial structure of the order 
parameter, we will now discuss how Major ana fermions 
emerge. An isolated vortex in the superconducting or- 
der parameter comprises a Majorana fermion as a bound 
state at precisely zero energy. However, in the vortex 
chain, these bound states can potentially couple to each 
other and form an energy band around zero energy. To 
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study this, we will construct an effective tight-binding 
model for the Majorana fermions in the vortex chain. 

First, let us consider an isolated vortex located at 
x = in the vicinity of a domain wall of type (ii), 
i.e., at \xi\ <C Z, in the strongly type-II limit, by as- 
suming <C Ao and \B\ <C Ao. We also assume 
that the vortex core radius, which is given by the co- 
herence length £ = v/Aq from the analogy to con- 
ventional type-II superconductors, is much shorter than 
the domain- wall thickness Z; £ <C I. Then, the prob- 
lem is reduced to a Bogoliubov-de Gennes Hamiltonian 
H A :=J d 2 x^{x)U A ^{x), where 



H vk A 

vk ] -[i A 

A* ji -vk 

\ A* -vtf ii 



(4) 



We used the bispinors Sfr^(x) = (^(x), — iijj T (x)a2) and 
the complex momentum operators k := —d\ + 182 and 
fct = d\ + 182- When B > 0, the order parameter A (a?) 
has the form A±(as) = ±A (|cc|) z/\x\ for |sc| <C £, where 
z = xi + and A (x) is a real non-negative function 
with A(0) = 0. With the sign ± we allow for two different 
choices of the gauge for A(x). A zero mode \£±(a?) bound 
to the vortex is the solution to T~L A± ^±{x) = 0, which 
is normalizable, i.e., J d 2 x^ ± (x)^±(x) < 00. To linear 
order in /i, it is given by 



1 

_/£_ 



\ 



v ±1 y 



exp 



1 Z 4 ^ 1 
^ Jo 



dxA(x) 



(5) 



where M is the normalization factor [19] . 

We turn to the case where A (a?) includes the vortex 
chain instead of an isolated vortex. The vortices are cen- 
tered at x = (n — \)x q , (n = 0, ±1, ±2, . . .). The overlap 
integral t of the Majorana modes bound to two neighbor- 
ing vortices can be written in terms of two solutions for 
the isolated vortex, and as 



t 



:= Jd*x*i(x+^)H A *-(x-^). 



(6) 



In our case of \B\ <€. Ao, which yields £ 7r/M> we can 
approximate A = iAo [A (a;) = Ao] in Eq. (rol. We obtain 



t = C[i 



Ao 
\B\ 



3/2 



exp 



ttAq 
2LBI 



O(fi 2 /A ) (7) 



as the leading asymptotic behavior in |/i|/A <C 1, where 
C is a dimensionless constant of order one that depends 
on the spatial structure of A(x). We note that the in- 
clusion of the in-plane component of the Zeeman field in 
Eq. Q for the calculation of t does not change the leading 
asymptotic behavior. The above example demonstrates 



that we can create a chain of Majorana fermions whose 
transfer integral t is controlled by /i as well as \B\/Aq. 

The vanishing of the hybridization and the emergence 
of many Majorana zero modes at \i = can also be 
explained in terms of Weinberg's index theorem on the 
number of fermion zero modes [20j [21]. Namely, there 
exist at least |n| zero-energy states for our Hamiltonian 
Eq. Q at fi — 0, where n := ^ § c dx-V arg A is the num- 
ber of vortices enclosed by a loop C. This property sur- 
vives, when the in-plane Zeeman field B(x) = M(x,0) 
is restored in the Hamiltonian Q, where M is in the 
form of Eq. (|3| with j = 3. Then, B at \x\\ —> 00 can 
be absorbed by a gauge transformation resulting in an 
exponentially decaying fictitious vector potential ^4, and 
hence the condition for the theorem [20] is satisfied. In 
contrast, if there remains a non-negligible out-of-plane 
Zeeman field B${x) unlike the above case, it couples the 
Majorana fermions and lifts their degeneracy. 

A tight-binding chain of Majorana fermions can be for- 
mulated in general with the real scalar Majorana fields 
4>j by the Hamiltonian, 



#eff = i + 5tj)(j)j(j)j 



(8) 



Here, j labels the vortices along the chain and Stj are 
random fluctuations of the hopping amplitude caused by 
spatial inhomogeneity of fi. A similar model was intro- 
duced by Kitaev in Ref.l22land has recently been studied 
in Ref. [23] Since the Majorana fermions are charge neu- 
tral, they do not contribute to the electric transport but 
to the thermal transport and the specific heat. The low- 
temperature one-dimensional thermal conductivity k due 
to the Majorana fermions is obtained within the semiclas- 
sical transport theory as n/T = 7t 2 /c!|£|t/(3M). Here &b 
is the Boltzmann constant and r the elastic scattering 
time associated with the backscattering due to the fluc- 
tuations 5tj. The presence of the Majorana fermions can 
be confirmed by the measurement of their T-linear contri- 
bution to ft, with the characteristic V-shaped dependence 
of n on I ji I across the Dirac nodal point. 

An alternative direct way of detecting the Majorana- 
fermion chain would be an STM experiment on the su- 
perconducting surface. This requires a careful fabrication 
of a spatial gap in the form of a strip separating two op- 
positely magnetized ferromagnets, to allow the STM tip 
to probe the superconducting surface. 

For these measurements, the temperature should be 
low enough to suppress contributions from other midgap 
states bound to the vortices. The second-lowest-energy 
midgap states appear at \E\\ = V^Ao/2 in the limit 
£\q\ —> and at /i = [24]. This leads to the condition 
k^T <C Ao. The suppression of phonon and magnon 
contributions to n also requires low temperature: The 
phonon contribution can be separated by its cubic tem- 
perature dependence. An insulating ferromagnetic thin 
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film with a large magnon gap can be used to suppress 
magnon contributions. 

We will now discuss a possible experimental realization 
of our proposal. If a 3D TI becomes superconducting 
upon cooling, the surface Dirac fermions associated with 
the 3D TI acquire an energy gap. Then, this may be de- 
scribed as an s- wave superconductivity of Dirac fermions 
hosting a single Dirac cone in the normal state. When 
an insulating ferromagnetic thin film is deposited on this 
surface, it will provide an ideal laboratory for our theo- 
retical model. Candidate materials for the 3D TI with 
superconducting instability include Cu^E^Ses which be- 
comes a strongly type-II superconductor by doping Cu 
into the parent 3D TI Bi2Se3 [25j[26], and the supercon- 
ductor LaPtBi which has been predicted theoretically as 
a 3D TI [27 j. The chemical potential needs to be care- 
fully controlled by applying a gate voltage or chemical 
doping. Insulating ferromagnetic thin films of spinel com- 
pounds CdCr 2 Se4 and CdCr 2 S4 [28] could be deposited 
on the surface. In particular, slightly off- stoichiometric 
compounds are useful for inducing appreciable magnetic 
anisotropy in the parent compounds [29] and thus sup- 
pressing the magnon contribution to the thermal con- 
ductance. A magnetic structure needs to be carefully ex- 
amined with Lorentz microscope. Typically, the domain 
structures are filamentary, thus creating several parallel 
chains of Majorana fermions. Then, their contributions 
to n may sum up to a larger signal. With realistic model 
parameters v « 0.3eVnm for the 3D TI Bi 2 Se3 [30 , and 
Ao ~ 0.5 meV, the spatial separation of two neighbor- 
ing vortices is given by n/q = ttv/2B « 5/im, provided 
B w 0.1 meV. The angle-resolved photoemission spec- 
troscopy measurement for Cua;Bi2Se3 at the optimal dop- 
ing x = 0.12 shows /i « 0.5 eV, which is about 0.25 eV 
above the bottom of the 3D conduction band [26 and 
clearly in the limit |/i| ^> Ao- In this case the transfer in- 
tegral t is expected to be much larger than the value esti- 
mated from Eq. (JT]) in the opposite limit, \t\ « 4xl0 _3 |/i|. 

Lastly, we propose an alternative way of creating a 
chain of Majorana fermions by using a heterostructure 
based on a Rashba- semiconductor, i.e., InAs, sandwiched 
by an s-wave superconductor and an insulating ferromag- 
net with a domain wall [Fig. [ljc)] . We apply a mag- 
netic field parallel to es to tilt the magnetization deep 
inside the two domains of the ferromagnetic insulator, 
yielding the effective Zeeman field in the semiconductor, 
B — ^ ±Bei + £?3e3, for x\ — ^ ±oo, where the domain 
wall extends along the e2-direction. The chemical poten- 
tial needs to be tuned closed to the T-point energy level, 
so that £3 gaps out the smaller of the two spin-orbit 
split Fermi surfaces. By linearizing the dispersion near 
the remaining Fermi surface, the effective noninteracting 
Hamiltonian takes the same form as Eq. ([!]). This can 
therefore produce a similar Majorana-fermion chain. In 
this case, however, one cannot largely control its energy 
bandwidth, since the chemical potential for the original 



electrons is always away from the Dirac nodal point. 

In conclusion, we have shown that a chain of Majorana 
fermions can be created along the domain wall of a fer- 
romagnetic insulator that is coupled to a superconductor 
with a single Fermi surface of helical electrons. In the 
search for devices that allow the manipulation of Majo- 
rana fermions as building blocks for a topological quan- 
tum computer, the proposed structure could in principle 
serve as a one-dimensional circuit path. In particular, 
the possibility to adjust the coupling between neighbor- 
ing Majorana fermions by changing the chemical poten- 
tial can be beneficial in this context. 
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